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The role of long range badly screened Coulomb interactions in a one-dimensional chain of Joseph- 
son junctions is studied. Correlation functions for the phase correlator are obtained as a function 
of the Josephson coupling energy, the short range part of Coulomb repulsion and its long range 
component. Though quasi-long range order is no longer possible and the usual Kosterlitz-Thouless 
, transition no longer exists, there are remnants of it. As an application, we calculate the I — V curves 

' for Andreev reflexion when a normal metal is placed in contact with the chain. Formally, there is 

always an offset voltage Vq below which no current can flow, however, in some regimes Vq can be 
' ^ , negligible. Contrary to what happens without long-range interactions, the Andreev current, as a 

' function of applied voltage, increases faster than any power law. Signatures of long range interac- 

' tions and phase slips appear in the I — V curves. Possible application for quasi one-dimensional 

thin superconducting wires is outlined. 
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I. INTRODUCTION 



U 

Oh; 

^ I Josephson junction arrays (JJA) have been widely studied, both theoretically and experimentally^. They show 
a number of phenomena, vortex physics^, quantum phase transitions. They are also a paradigm; a large variety of 
' physical systems can be modelled using Josephson junctions arrays, for example, one-dimensional thin superconducting 
wires^ and granular superconductors. Nowadays, they are used to realize protected qbits* or for metrological 
' applications, for establishing current standards through frequency measurements of Shapiro steps^. 
''O . Here, we shall restrict ourselves to the 1-dimensional case. Furthermore, we shall neglect dissipation, which is 
Ch ■ crucial in some quantum wires. No coupling to a dissipative environment is considered. For normal regular junction 
chains, the essential physics is the competition between Josephson coupling and local Coulomb interactions. In one 
dimension, true long range order can not exist, the correlations of the superconducting order parameter decay at best as 
a power law. When local Coulomb interactions are increased, the power law is modified, and eventually a Berezinskii- 
Kosterlitz-Thouless (BKT) transition occurs. Beyond the transition, when Coulomb interaction dominates, the decay 
of correlations becomes exponential. This is due to the unbinding of phase slips, whose existence have been shown 
directly in recent experiments^. 
' Topological defects have always played a central role, from BKT transitions in magnets, roughening of interfaces, 
, disordered Luttinger liquids and some stripe models for high-Tc superconductors. 

■ We would like in this paper to study the role of long range Coulomb interactions, hereafter denoted as LRI, which 
' have not been taken into account so far in one-dimensional JJ chains. One motivation is that, in one-dimensional 
i: clean superconducting wires, the long range part of the Coulomb interactions is poorly screened. For example, the 

T—{ ' plasma mode frequency goes to zero for small frequency-. This has motivated some recent studies^'. Clearly, in these 
L| i; systems, fluctuations of the phase of the order parameter have to be taken into account. The Hamiltonian looks very 
. 1-H ■ much like the one of a one-dimensional Josephson array, with a long range Coulomb term in addition. 

, Long range interactions, hereafter denoted by LRI, were studied a long time ago at the spin- wave level for electronic 

■ systems2.ii£. They also appear naturally in modelling underdoped high- Tc superconductors, where phase fluctuations 
of the order parameter are importan t^^'^^ . Influence of LRI have also been studied recentlyi^ theoretically in two and 
three dimensions. Our goal is to obtain the phase correlator of the superconducting order parameter as a function of 
time and space. This quantity governs numerous physical properties. In this paper, we shall give only one application. 
Namely, we calculate the average Andreev current from a normal metal to a JJ chain. Because LRI are normally 
screened in fabricated JJ chains, the model discussed here is not meant to be a realistic model for a JJ chain but 
can be relevant to thin 1-dimensional superconducting wires. However, calculations developped in this paper can be 
adapted to take into account the screening of the Coulomb interaction between islands and make predictions about 
Andreev current in real JJ chains. 

After the introduction, the body of the paper is organized in four parts. 

In part two, we first describe the model of JJ chain we use, and its connection to high- Tc materials. It is convenient 
to go to imaginary time. When the Josephson coupling dominates over the Coulomb repulsion, if phase slips are 
neglected, the modes are spin waves. The long range interactions cause the decay of correlations to change from those 
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of a Luttinger liquid to a more pronounced decaj*is. 

In the next section, we take into account phase slips. The Hamiltonian decouples into a part describing spin-waves 
(long wavelength excitations) and another part describing phase slips (rapid variations of the phase). In the last 
sections, we discuss the interaction between phase slips. For distances much smaller than a characteristic distance, 
they interact logarithmically. Beyond this distance, the interaction behaves only as the square root of the logarithm 
of the distancei^. 

In part three, we examine how the KT real space renormalization group (RG) analysis is modified by LRI. There, 
as soon as LRI are present, strictly speaking, no quasi long range order is possible anymore. However, the initial long 
range interactions can be very weak. If we iterate the RG flow, it will take some time for LRI to alter the behavior of 
the usual KT flow. The LRI will show up only after a certain crossover "time" , which we calculate and corresponds 
to some temperature T* . On the other hand, RG has to be stopped at some cutoff "time", H/kBT. Therefore, in 
order to see the influence of LRI, the temperature must be lower than T* . The dependence of T* on the values of 
parameters such as the Coulomb energy, Josephson coupling and strength of LRI is explicited. 

Part four is devoted to the calculation of the phase correlation function. Analytical continuation is required to 
obtain the correlator in real time. This part is rather technical; some aspects of the calculations are relegated to 
appendices. The main result is the expression of the phase correlator, which differs from the usual Ornstein-Zernike 
expression. 

In the last part, we discuss one physical situation. The Andreev current between a metal and a JJ chain, as a 
function of applied voltage is calculated. We compare the curves obtained in three cases. The first case corresponds 
to only on-site Coulomb interaction. In the second case, a Coulomb interaction between nearest-neighbor islands 
is added. This model is often used in the literature. Finally, in the third case, LRI are added. For relatively low 
voltages, the Andreev current remains always smaller with LRI. We focuse on the case where, without LRI, we would 
be in the disordered phase. There is a threshold voltage Vq below which no Andreev current can fiow. This remains 
true in the presence of LRI. However, the current versus voltage (/ — V) curves are different from the case without 
LRI. Without LRI, the Andreev current starts as {V — Vo)"^^^, whereas with LRI, it starts to increase faster than any 
power law. This can be viewed as a signature of LRI and could be tested in experiments. 

In the conclusion, we discuss two possible applications, one for realistic fabricated JJ chains and the other for 
one-dimensional superconducting wires. For JJ chains, the interaction between island charges is screened but the 
methods used in this paper can be adapted to make at least qualitative predictions. For Id wires, owing to the strong 
dependence of a phase slip fugacity with temperature and wire width, and also for other reasons detailed below, it is 
not obvious that LRI can be detected experimentally. Maybe they could for very long wires. 

II. GENERALIZED VILLAIN MODEL AND EFFECTIVE ACTION WITH LONG RANGE, XY MODEL 

IN TWO DIMENSIONS 

A. Imaginary time action 

We would like to study the influence of poorly screened Coulomb interactions in 1 dimension. One possibility, which 
has been widely used some time ago, mainly in dimensions two and three, would be to take the Hamiltonian for a BCS 
superconductor— li^i and to introduce Hubbard-Stratonovich transformations in the particle-hole and particle-particle 
channels in order to make the usual gap A — |A|e"^ appear. The amplitude fluctuations are integrated out, so are 
the phase of the density fluctuations, leaving us with a phase action only. This procedure has been used in Ref. |^, 
except that they stayed at the RPA level and did not consider topological excitations for the phase (p of the order 
parameter. Another way would be to start with a one-dimensional array of Josephson junctions, modelled by an ideal 
Josephson junction (JJ) and a capacity Cj in parallel. In addition, on each island, the capacity to the ground has Cg 
to be taken into account. See Fig. 1. 

The Hamiltonian is 

HjjA = XI ^ cos(6l,+i + ^i^id^j 

with qi = 2eni and the interaction potential Ui,j — 2e^C^j and C^j is the capacitance matrix. The last term of Eq. 
(1) represents LRI, A is a constant and C a capacity. Usually, only the capacity of one junction Cj and the capacity 
from an island to the ground Cg are taken into account, the capacity matrix Cij has the form Ci,j = Cg -\- 2Cj if 
i = j, Ci,j = —Cj if i = j ± 1, and Cij = otherwise. It is reasonable to approximate the interaction potential^ by 
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FIG. 1: One-dimensional array of Josephson junctions with long-range Coulomb interactions. 



the form 



2J 



c. 



(2) 



with A — ^Cj jCg. The screening length is of the order aA, with a the interjunction spacing. Typically, for one- 
dimensional fabricated chains^! Cj ^ Cg but still A is of order 10 to 30, not more. LRI thus do not appear to be 
completely realistic for fabricated JJ chains, but for thin superconducting wires, the Coulomb interaction is badly 
screenedi2ii2^ and can be longer range. We shall thus take a Hamiltonian of the form 



HjJA = ^^J cos{0. 
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(3) 



with C ^ typically an effective capacitance to the ground. We can also see this Hamiltonian as a Id analog of the 
following Hamiltonian for superconductors with fluctuations of the phase of the order parameter 



^ 2m* 



E 



1 m*v'' 
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(4) 



In this case, is the superfluid density in 3d. p is the charge density, 9 will denote the phase of the order parameter. 
Its conjugate variable is H, which is proportional to the charge q. When phase slips are neglected, the procedure to 
obtain an effective action for 9 is standard and we just recall it here for completeness. 
Writing qi = 2eni with rii integers, 9i and Ui are conjugate variables, [ni,9i] — Sij, 



H = ^-Jcos{9, - 



(5) 



with = ^ the charging energy for nearest neighbor. Following Ref. [T^ T denoting the temperature and 
/3 = l/{kBT), the partition function is calculated in imaginary time, by discretizing in time steps of size H/ ^/2EoJ . 
The system without LRI would be equivalent to a usual XY model with isotropic coupling. In the following, 9ij will 
no longer be operators but numbers; index i from 1 to A'^ = L/a denotes the position in real space and the second 
index j from 1 to iV' = I3\/2EqJ the position on the imaginary time axis. Making h = 1, the partition function reads 



Z = ^exp 



2En 



cos 



cos 



(6) 
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where the sum is on all configurations {9}. 

Now, in the spin-wave approximation, the cosines are expanded and the Lagrangean L and action S take the form 

L dr. (8) 



Applying a Fourier transform and integrating out the Ilq ~ Uq field. 



1 



4.Eq + AAEoI{q) 



dr, (9) 



where I{q) is the Fourier transform of the Coulomb potential, regularized^ at very small q. dr9q is in fact — 9q j , 

and 9qj is the Fourier transform in space of 9i j . 

Rescaling the imaginary time and counting it in units of h/\/2EoJ, and setting h — I gives the effective action on 
a discrete space-time lattice 

^ ^' 1 I T~ 2 I 7~ 1 ^' 



with N' — (3^/2EoJ and 9 the Fourier transform in space of 9ij defined as 

N 

1 

^9 ^ 



1 ^ 

= ^ E (11) 



A^ 1 
n— 1 

for g multiples of 2tt/N. 

Now, taking into account the periodic character of 0, we can take the effective action 

N N' I— I- N- 

Seff = E E V 2^ - ^-^■) + V 8^ ^ YTAIiq) T.i.^"^^^^ - ^.^.) (^-.-.+1 - ^-.^.)- (12) 

Having derived an effective action, we first turn to the case where the periodic character of the phase can be 
neglected. The approximation is good for large J/Eq. Cooper pairs feel very little Coulomb-blockade. The fugacity of 
phase slips goes essentially^ as exp(— 4-^/2 J/i^o). Phase slips are neglected and a phase correlator derived. Owing to 
LRI, the system is no more a Luttinger liquid. Correlations ressemble some correlators appearing in one-dimensional 
quantum wires with low electron densities^iiiS, sometimes dubbed "Wigner crystals" . 



B. Spin waves 



In this section, we allow only smooth variations oi 9. It is then possible to expand the cosine term. This gives 
immediately a spin-wave Hamiltonian. We check here the known spin wave dispersion relations. 



5, 



SEq 



1 



1 + AI{q) 



(13) 



with I{q) ~ ln(g ^) for small q. Setting K ~ ^j"-, this results in 



^q,uj^ — q,—uj I 



1 + 



and thus 



G{x,t) = {{9{x,t)-9{0,0))') = 2^J^-E / 



1 — cos(cj„T -f- qx) dq 



2 I 

1 l+AI{q) 



(14) 



(15) 
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where /3 = \/{kBT) and w„ are the bosonic Matsubara frequencies. For distances x and times r not too long, so 
that the long range term does not influence the behavior, the correlations are those of a Luttinger liquid G(a;, r) — 

K In {^^^-^^^—^ 1 where a is a short range cutoff (corresponding to finite bandwidth), whereas for very large timesi^, 



G{x,t) = kV2A In^/^ 



if T is counted in units of h/y/2EoJ. The correlation function of e*^ for spin-waves read 



/2ATVlnr 



je{0.0)-i9{Q.T)\ _ 



exp —B 



ln(T) + 2 ^'^(l^ 



3/2^ 



(16) 



(17) 



with B = Ky We would like to calculate the full correlation function, not only in the spin wave approximation, 

but including phase-slips. This turns out to be more intricate, we will neglect all factors that vary slower than any 
power law, such as, for example, logarithmic corrections or terms as exp(— CoVln r) with Co a constant. Within this 
approximation, we have the spin-wave correlator 



g(r) ~ exp 



-B 



Mr) 



-1 3/2N 



(18) 



C. Phase slip Hamiltonian 

Now, we take into account phase slips, that is vortices for the XY model in 1 space and 1 time dimension. The 
procedure to split the action into two parts, one due to spin- waves and the other to vortices is standard without LRI, 
however we recall it here with some details to point out the influence of LRI. 

We use the model of Villain and transform the exponential of the cosine term in as the sum over n integers 
but with a quadratic actionii, involving and n. The space direction will be denoted as x and the imaginary time 
direction as y. The long range interaction occurs only in the x direction, for bonds directed in the y direction. Let us 
denote by 



As a consequence, the Coulomb term, i.e. the second term on the r.h.s of Eq. (|T2)) . has to be modified from 
E,^^(9)Ef=i(^.,.+i-^,,.)(^-...+i-^-.,.) toj:,F{q) j:,[Ei{&U+i-Oi,j-nl)e^'^'] [E.(^r,,+i -0.,, -<,)e''^'-^ , 
where denotes the (27r times) integer attached to the bond which goes from site (r, j) to site (r, j+1), and pertaining 
to the y direction. The ^ are simply the Fourier transforms of the njf,„. The lattice spacing will be denoted a, 
and e^^ (e^) will be the unit vector in the x (y) direction. 

The part not involving the n^'s is of the type J2q -^il)J2k ^q-ky 0-q,-ky so that finally the total action is of the 
form 

J2 I^^^a0k-27rn£|2 + e|^^fe|^ (20) 

a—x,y k 

where 

= 2sm(k.a"/2), (21) 

and = ae^, — aey where e^; and Gj, are the unit lattice vectors, e is an infinitesimally small positive quantity. 
From now on, we will denote q = kx so that k has components {kx,ky). The only difference with the normal XY 
model is that the coupling constant in the y direction has a dependence in k^ for small k^ (long wavelengths). 

Setting Ax ~ J 2^ a-iid -^y = -^xFiq), 



k a—x^y 



(22) 
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Diagonalizing this action yields the decomposition into a continuous field 6k and a discrete field rik- The action must 
be of the form 



Seff = Se + Sv, 

k 



(23) 
(24) 

(25) 



where a and 7 take two values, x and y. The quantities Da^-f{k) and C(fc) arc functions of k. The linear transformation 
is given explicitly by 



Ok + 2i7r 



Da,j(k) = 47r 



a—x,y 
2 



-^K-?^^ + -^yKy 



KaK^ (2(5a,7 ~ 1)' 



(26) 
(27) 

(28) 



and Sa^-y = 1, if Of = 7 and zero otherwise. 

Next, Sv is rewritten in terms of vortices defined as 



which yields the vortex Hamiltonian 



fcT^O 



qkQ-k 



(29) 



(30) 



The apparent divergence at fc = is in fact absent because of charge neutrality {qk=o — J2i 1i — 0)- Going back to 
real space, denoting Ntot = N N', the total number of sites: 



V J^tot y. 



enables to rewrite the vortex Hamiltonian as 



(31) 



(32) 



with, taking into account values of Ax, Ay and F{q), 



P,P' 



[l + AI{k^)]Kl+Kl 



(33) 



First, the self-interaction V^.p has to be withdrawn, but it is not so crucial because of charge neutrality (X]p = 0)- 
Second, the vortex Hamiltonian as it stands is incomplete. The vortex fugacity term, of the form ^pi^a y)qp is missing, 
with In y being the fugacity of a vortex. In fact, one should start with a generalized Villain modeli^, introduced by 
Jose, Kadanoff, Kirkpatrick and Nelson, hereafter denoted by JKKN. As noted in Ref. [l^ such a term is naturally 
generated by real space (Migdal-Kadanoff) renormalization. It is therefore necessary to reinstall this term. Next, we 
look at the shape of the interaction between vortices. 



D. Two regimes 



Although it is known that, with LRI, the vortex-vortex coupling for large separations^ (much larger than the lattice 
constant) goes only as in -^/In |Rp — Rp'|, we need to derive a quantitative expression. The behavior is dominated by 
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the behavior of the integral for smaU |k|. For theses small wave- vectors, Kx — ak^ and Ky ~ aky, a being the lattice 
constant. We can take Rp — and Rp/ = p cos(z/)e2. + p sm{i')ey. 



I = 



J 
2^ 



1 — cos 



(jtxp' cos(i') + kyp' sin(j/) 



1 + A/(fc^))/c2 + 



(34) 



The cosine term will be close to 1 only for k^ much lower than l/[p' cos{iy)], and ky much lower than sin(^)]. To 
extract the main behavior for large distances p', we can replace it by the integral 



J 
2E'q 



1 



(p' cos(i/)) (p' sin(i^)) -'- 



{1 + Aln(fc^i)}fc2 + fc2 



(35) 



where km is a wave- vector much larger than y , but still sufficiently small that the expansion of and Ky for small 
k holds. The integral reads 



J 



at an 



atan 



^A:^y^l + Aln(/s:r^)^ ^p'sm{l^)kx^J^^^Ah^(J^)^ . 



(36) 



As p' gets large, for different from and 7r/2, p' cos{i/) / [A^/]n{f7)\ goes to infinity, so the first atan term in the 
bracket is close to tt/2 whereas the second atan term remains much smaller than 1, provided tan(i^) is not too small. 
We can thus safely replace the term into brackets by 7r/2. The integral on k^ is then readily performed and is in 



'l + A In(fc-i). The limit A = is singular. It is better to discuss limiting cases. 
Case 1: • ln(^p' cos{v)j < 

This corresponds to distances where the decrease of the correlations due to long range unscreened Coulomb inter- 
actions still has not taken place. We have then 



p,p' 



2tt 



cos(i/) 



(37) 



Case 2: • In^p' cos(i^)^ > 

This corresponds to distances where the decrease of the correlations due to long range unscreened Coulomb inter- 
actions has taken place. 
We have then 



27r 



^yin(p'cos(.))M. 



(38) 



Thus, at very large distances, there is no longer sufficient interaction between vortices to hold them tied together, 
whatever the short range term is^^. Mathematically the BKT transition is expected to disappear. However, we need 
to put this on a more quantitative basis. 

Furthermore, in real systems, the Coulomb interactions remain unscreened only up to a certain screening length 
^ej which can be large but is not infinite. Also, the magnitude of the LRI Coulomb term may be much smaller than 
the one of the local Coulomb term. In order to determine the phase diagram and the phase correlations, we need 
to implement a real space renormalization group scheme analog to the one developped in Refs. [2IIT8I. Care has to 
be taken since LRI will eventually dominate under the RG group. The last term in Eq. (2) is more relevant than 
the second one. However, under the RG fow, it will take some time for the last term (LRI) to overcome the effect 
of the local one (i.e. the second term in Eq. (2)). The situation is reminiscent of LRI in one-dimensional quantum 
wires with very low electron densities; there are important differences however. There will be remnants of the BKT 
transition. To proceed, we perform a RG in real space closely following Kostcrlitz'original paper—. 
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III. RENORMALIZATION GROUP IN REAL SPACE 



A quick look at the situation is that, owing to the Coulomb potential, at very large distances R, the interaction 
between vortices goes like y/Tap (instead of In p in the usual case) , and is too weak to overcome the entropy term 
for placing a free vortex, which goes as In p. Nevertheless, the last terms in Eqs. (|3I4I) can be small compared to the 
second ones, even if they are less relevant in RG sense, their initial value can be much smaller. Thus, it makes sense 
to treat them as a small perturbation, at least initially. This leads us to investigate how the usual KT RG flow is 
modified. 



A. Modified KT equations 



Here, the partition function has a form which is slightly different from the usual KT flow. The gas is neutral, 
without electric field, n will be the number of "+" charges. 



^ = Ei^'" exp{^^^(r,-r,),,q,ln|^l-^| } 



where J' is a coupling constant, — 27r^/ and is a function that behaves as 1 for not too large r,; 
decays as In jr^ — rj]^^/"^ for very large 



2Et 

Ti — rj\. A model T can be 



F{r, - r^) 



1 



A In Wj - 



(39) 



but 



(40) 



T is the short time cutoff and K, the fugacity. 

In the case where the long range term is much smaller than the short range one (at the beginning of the renormal- 
ization procedure only), we can however follow the same steps as Kosterlitz. The details are shown in Appendix A. 
We express the RG equation, using the variables 



X 
Y 

z 

Now, with the long range interaction, this becomes 

dX 

dz 
dY 

dz 

with B ~ 

a line of (stable or unstable) fixed points. 



In r. 



-Y - B^z, 
-2XY, 



(41) 
(42) 
(43) 



(44) 
(45) 



A. Note that B can be much smaller than one. Now, because B is non-zero, the line F = is no longer 



B. Physical description of the flow 

The flow is no longer autonomous, as soon as S ^ 0. We expect a deformation of the usual flow. In particular, the 
line Y = which was a line of fixed points in regular KT flow is no longer a true line of fixed points. 

If the LRI were absent, we could be either in the case where phase correlations are decaying exponentially (disordered 
phase) or in the "ordered" phase, where correlations decay as a power law. In the first case, adding the long range 
interactions will just make the phase correlation length smaller. The more interesting case is when we start from the 
ordered phase. We thus examine quantitatively how switching LRI drives us into the disordered phase and, in the 
next section, we shall examine what the phase correlations look like. 

In the usual case, (without LRI), the flow is autonomous (i.e. does not depend explicitly on time) and trajectories 
are of the form X^ — Y = Xq —Yq = C, where (Xq, Yq) is the starting point. Setting y = VY, they are hyperbolae 
in the {X,y) plane. In the "ordered" phase, the flow on the invariant curves is: 
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FIG. 2: Modified KT flow in the presence of long-range interactions. Only a point which would correspond to the "ordered" 
phase of the usual KT diagram is shown. Here, because of long-range interactions, we always end in a completely disordered 
phase. 



which integrates to 



2 




In 



2(z-zo). 



(47) 



zq is the initial "time" which can be set to zero. Integration to obtain X is also straightforward. 

The important thing is that, as z becomes very large, Y , which is related to the fugacity, tends to zero exponentially, 
while X tends to a strictly positive value X/i„, see Fig. 2. 

In our case, (with LRI), we can see what happens by performing the following approximation: 

- if By/z is much smaller than Y , we neglect it, 

- when, it becomes of the order of Y , we do the opposite and neglect Y in the r. h. s. of the dX/dz expression, 
Eq. dllD. 

Thus, there are at least two stages. In the first stage, we have the usual KT flow whereas in the second stage the 
long range term comes into play. We shall see that there are three stages actually. We now describe the first stage. 
Starting from a definite value of Y , by the regular KT flow, Y will decrease. Suppose that Y has decreased by a 
substantial amount, such that Y will be much smaller than C. This enables to replace ^1 -|- y^y by ^ ylz) • Inverting 
the relation yields 



Y{z) 



4C 



exp(-VC{z - za) 



(48) 



with D = JS- ' ' ' 



Yo ' V "'^ ^ necessarily much larger than one). The physics is that C is proportional to 

the initial distance to the separatrix. It is proportional to (Tc — T)/Tc m the thermal X-Y model. Now, we have to 
stop when Y reaches the value B^Jz — zq, because we are no longer allowed to neglect the term A^/z in Eq. p4)) . 
This happens when 



4C 

zJ 



exp 



-VC{z-zo)) =6/ 



Zo- 



(49) 



To get the main feeling, we can show that the typical value of z for which this happens is zi , such that zi — zq — ■ 
This corresponds to a typical imaginary time ri, or corresponding lengthscale li such that 



n 



h = 



Neglecting In(ln) contributions, Ya, the corresponding value of Y is 



Yr, 



(50) 



(51) 
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Ya goes to zero as B goes to zero. Now, we will enter the second stage. 

We now describe this second stage, where the second term in Eq. becomes larger than Y. Then, for z larger 
than Za, we integrate the RG equations, neglecting the Y term with respect to B^/z. 

X = Xa-^z'/', (52) 
and, using the other flow equation to obtain y, 

Q 

\nY = \nYa-2XaZ + —Bz^/'^. (53) 

Then, X will decrease and Y will continue to decrease, till X passes the value 0. Y will be minimum ioi X = 
and then will start to increase. In fact, it increases faster than exponentially, so it will quickly become larger than 
[By^'z-^^'\ , so that our approximation will fail. The value of Z2 for which this happens can be estimated, but anyway 
Z2 — zq <^ zi — zq, corresponding to a length I2 = e^^ ^ h. 

Then, there is a third stage, where Y becomes larger than [\/B{z — zq)] • Then, we can again neglect [^/B{z — zq)] 
in front of Y. We are then in the third regime which has analogies with the KT flow but on the disordered part of the 

transition. This lasts till Y has reached an appreciable value such that one is no longer in the vicinity of the critical 

i_ 

point. Let Z3 be the corresponding value of z, and I3 = e^^~^^ . I3 will be also of order B . Since C is proportional 
to Tc — T, where Tc is the transition temperature without LRI, it is tempting to conclude that, since B goes as 
the correlation length with LRI f ' behaves as 

h 

with D an unimportant constant of order one. It has the same form as the KT correlation length in the disordered 
phase (without LRI) but \n{A~^) replaces the usual^ constant b. Of course, tends to infinity as B tends to zero. 
The general shape of the "flow diagram" is best seen on Fig. 2. 

Renormahzation has to be stopped when r becomes equal to ph. Equating ^' to h/{kBT) yields a temperature T* 
given by 

\n{A- 



e = D^^e^p(\n{A-')/2VT~f), (54) 



T* = iDkB)-'V2Ej e^p\^-^^^=^j, (55) 

with Tc the critical temperature in the absence of LRI. If experiments are performed at a temperature larger than 
T* , it will not be possible to see the influence of LRI. 

Having the general shape of the "flow-diagram" , we now turn to the determination of the phase correlation function. 



IV. CORRELATION FUNCTIONS 



This section is a little technical. Just as the Hamiltonian decouples into a vortex part, the phase correlation function 
also do. Estimation of the spin-wave part is straightforward, but the estimation of the vortex part is more intricate. 
In section A, we first determine the contribution of one vortex to the phase correlator. This involves a vortex influence 
function v{r). Comparison between the usual case and the case with LRI is emphasized. 

Then, in section B, average must be carried out over the possible conflgurations of vortices, using the Boltzmann 
factor with the vortex Hamiltonian. A Coulomb gas formulation is necessary. Distinction must be done between 
smooth variations of the vortex influence function v{r), studied in section C, and abrupt variations (section D); there 
are thus two contributions. Summing both contributions gives the vortex part of the phase corrrelator. Adding the 
spin-wave part yields the total phase correlator in imaginary time. In section F, analytical continuation is performed 
to obtain it in real time. The phase correlator in real time is necessary for the calculation of the transport properties 
of J J chains. 

For the determination of the correlation function, we use the method of JKKN. In the vortex Hamitonian, the core 
energy has to added. Such a term naturally arises from their Migdal-Kadanoff renormahzation group. Therefore, we 
must use the generalized Villain lodel. The corresponding vortex Hamiltonian is 

Hv =Y. ^py'ip 1p' + (In J^) E 9P' (56) 
p.p' p 

with y the charge fugacity. The qp are integer charges living on the dual lattice of the original lattice. They represent 
vortices. The correlation function to be computed is (e'®*^^''^^e~*''^°'°''). Written in Coulomb gas language, this 
correlation function disentangles into a spin- wave part and a vortex part. 



11 



A. Determination of the vortex influence function 



We closely follow the methodi^'^S of Pelcovits. The first thing to do is to determine what is the influence of a vortex 
sitting at R on the vortex part of the correlation function .gy(r — r'). The total correlation function is 



gtot{r 



(57) 



Following JKKN, gtot{^ — r') is written as a product of nearest neighbor correlation functions, along a path V that 
connects r to r' on the real lattice. This is for the total correlation function (vortices plus spin- waves). After 
(Gaussian) integration of the spin-wave degrees of freedom, 



gtot{r 

The vortex part, gv{^ ~~ r') is given by 



') = gsw{r-r')gv{r-r'). 



with 



gv{r-r') = (^exp(^^m(R)^;(R))^, 

R, 



v{R) = ^V^(R-R')[^'(R)-^'-(R')], 



(58) 



(59) 



(60) 



R' 



where R are all the positions on the dual lattice, the function ?7'(R) {t]^(R)) is zero everywhere except on sites of the 
dual lattice that lie immediately on the left (right) of the path chosen joining r to r'. V(R — R') is simply Vp^p', see 
Eq. (p4)) . The brackets mean thermodynamical average with respect to Hy- 

This results directly from Gaussian integration, it can be understood as follows: the contribution of a vortex sitting 

at location R, and vorticity m(R) is exp^i m(R) w(R)^ . The function u(R) depends on r and r'. It has also another 

physical interpretation. It is the phase generated at site R by a vortex- antivortex pair, where one vortex has m ~ 1 
and is sitting on r and the antivortex (m = — 1) is sitting on r'. 

The determination of w(R) has been done in the literature without the long range Coulomb interaction (LRI). This 
function will be called w(R), and the potential energy between vortices will be denoted as V^p,. 

We explain how the LRI affect the results. Following JKKN, without LRI, 



(61) 



where G"(R-R') = j^, 



R' 



V°p. See Fig. 3. 




ri^ (R') - ri''(R')^ here 



FIG. 3: Determination of influence of vortex sitting at R on the correlation function (e* ^"''e 



ie(r)„~ie(r')\ 



For not too large distances, the interaction between vortices is still logarithmic and we just recall the known results 
for completeness. For simplicity, we take a path joining r and r' to be a straight line (see Fig. 3). Working in cartesian 
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coordinates, and taking the origin at r, the x axis, along the line going from r to r', the coordinates of R, r and r' 
will be {x, y), (0, 0) and (r', 0) respectively. The current integration variable R' on the dual lattice immediately above 
the cut will be of coordinate (x', y'). 



(R) = liniy'^o I -j-f In^^i/ {x - x'Y + {y - y')^) dx' = atan(2//a;) - atan(?//(a; - r')) , 



dy 



The function m(R) at large distances must be replaced now by t;(R). 



2tt jjr r d 

Integrating by parts, to extract the main contribution at large distances, yields the result 



(62) 



which gives 

u(R) = $(r - R) - $(r' - R), (63) 

with $(Z) = dXmi{y^lx,^), if Z = {xz,yz)- 
Now, consider the LRI, at large distances 



(64) 



«(R) = Hruy,^,^^!^^ J^^ ± (^^ln(y(^^2 + (y_y,)2)^ ^x' . (65) 



= — 7^\T^\ ; ^=atan( ) , atan(-) 

^-1/2 r-p 

' F, (66) 



V2 V 2£;o 

with F a quantity which is smaller by a factor ^ln(a;^ + y^)^ than the main one. Thus, v(R,) behaves^ as ln(R) 
for large R. 



B. High temperature expansion of the vortex part of the phase correlator 



From Eq. the average of exp(^i m(R)i;(R) j has to be performed with the full vortex Hamiltonian (in- 

cluding core energy). The Boltzmann weight of a particular configuration m(R) is given by exp^^(m(R))^ with 
A{m{R)) = Er^'W + Er,r' l^R-R' and thus 



, E{™(R)} exp A m(R) + z Er m(R)T;(R) 

9v{r-r')^ ^ ^- ^ . (67) 

E{m(R)}Cxp(^A(m(R))j 

A high temperature expansion of gy(r — r') is needed. Applying the Poisson formula fifn) = 
X^^-oo /('?^) ^^^"'"^ ^^'^ function /(m) = exp^A(m(R)) + i Er ™(R-)^'(R-)) and integrating over the 

continuous auxiliary field (j) enables to write it as 

gvir^r')^^, (68) 

with 

Z' = ""M^l + 2^pW] H'^') + 2'rp(R')] (TO(0)m(R))co„*') , (69) 

{p(R)} ^ R,R' ^ 
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where the sum is over all configurations of p(R). {m{0)m{'R)) cont is the correlation function calculated with weight 
A(m(R)) but as though the m were continous variables. Z is the same as Z' except that w(R) is absent. However, 
because a switch from the m(R) variables to the p(R) variables (through the Poisson formula) has been made, now a 
low temperature expansion is needed. A high T expansion in the m(R) variables transcribes into a low T expansion 
in the p(R) variables. In the absence of LRI, (m(0)m(R))cont has been calculated in Refs. [2lll2^ Even with the LRI, 
the correlation (m(0)m(R))cont remains short range and can be approximated by nearest-neighbor only. Substituting 
the exact (TO(0)m(R))cont by a nearest neighbor interaction results in 

Z' ^ J2 ^^P(~l^\/^ ^ {viR)+2np{R)-u{R')-2np{R'))'). (70) 

{p(R)} ^ (R,R') 

Following Refs. [T9ll20l a low temperature expansion is performed, keeping only the essential diagrams. To simplify 
further, only the non-overhanging self-avoiding-walks (SAW's) joining R to r' are kept which results in 

gv{r-r') = ^ P„(r - rV"/^^ exp(^^y^ J f \W nv{R)\' d'R^ , (71) 

where J J' means integration over all the plane except on the SAW. -Pn(r ^ r') is the number of non-overhanging 
SAW's of n steps joining r and r'. 



C. Smooth variations of w(R) 

The calculation of the correlation function splits into the integration on the smooth variations of the field v(R) 
and on the abrupt discontinuities caused by the different diagrams occuring in the high T expansion of Z' . In our 
case, the situation looks a little more complicated because v{R) does not satisfy the Laplace equation. However, this 
causes no trouble because the volume contributions coming from Aw, where A is the Laplacian, will turn out to be 
very small. The main contributions come from the regions near the path V going from r to r'. 

More precisely, the weight 



W = exp 



1 

8^ 



2Eo 
J 



J I" (Vr«(R) 



d^R 



(72) 



is to be evaluated, where the prime means integrations everywhere in the plane except on the diagrams. 
In the usual case, integration is made by parts 



(vRM(R))'d2jj^ = T1+T2 + T3, 



Ti 



u 



Au d^R, 



s 

u Vu .da. 



(73) 
(74) 
(75) 
(76) 



Tin 



S denotes the shaded aera in Fig. 4, Tin is a contour which encloses the cut (from r to r') at a small distance 60 and 
Tout is a circle of radius TZ where TZ will eventually be allowed to go to 00. d<j is the unit normal to the graph. On 
Fig. 4, it is simply the normal to the line joining r to r'. The first term Ti vanishes since Am = 0. The second term 
T2 will tend to zero as TZ goes to infinity. The last term has to be calculated using Eq. (pT|) . Care has to be taken 
that G"(R — R') behaves only as ln(R — R') for sufficiently large distances and this form of G' should not be used 
for small distances. 

With this in mind, taking r = and thus r' = |r' — r|. 



atan( — ) — atan( — ^ — 
x' r — x' 



dy' 



atan(-^) — atan( — 



y 



dx . 



(77) 



The factor 2 comes from the fact that integration above and below the cut r,r' has to be performed (see Fig. 4). 
We take the integration path a finite amount 60 above (or below) the cut to avoid spurious divergencies (see above). 
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FIG. 4: Determination of influence of the smooth variations of the field «(R) on the vortex part of the correlation function 



Then 60 will be unimportant. The angle difference in the difference of the two atan functions Eq. (177)) will correspond 

to the change in angle, taken to be tt. Dropping the spurious ln(6g) terms yields^ exactly T3 = 47rln|^|r' — rj^ 

Now, we turn to the case of LRI. The analog of the terms T^, i = 1 to 3 will be denoted as T/ and u has to be 
replaced by v. Calculations are tedious and relegated to the appendices. It turns out that T{ and T2 are completely 



negligible whereas Tg = Bi In 



Inlr 



, with Bi a constant. Even after taking the exponential of this will only lead 



to logarithmic factors for the correlation functions and we thus omit this factor. 



D. Abrupt variations of w(R), graph contributions 



Here, the part due to abrupt variations of the field w(R), is the same as without LRI sum on SAW joining r and r' 
and the result is analogous to the usual case, namely, the contribution, called gv saw reads: 

exp(-^^) 

gvsAw - /, ,1 (78) 

ujp is the plasma frequency ^/8Eq J /H. 

Being the result of a high-T expansion, this expression corresponds to short correlation length ^. It is probably not 
valid for very large In the limit A — 0, the usual divergence close to the BKT transition must be retrieved. 



E. Final form 



The spin- wave part, the vortex part coming from smooth variations of w(R) and the vortex part coming from abrupt 
variations of w(R) just need to be multiplied. Thus we have, with r — r' = |r — r'|. 



g{r - r') = cxp(^-B[ln\{r - r')ujp]\^^^^ exp{Bi In \n\{r - r')ujp\) 



exp(- 



|(r-r')^p 



(79) 



^(r - r')Ljp 

with Bi a constant. Contrary to the usual case, the part of the vortex contribution corresponding to smooth variations 
of v{R) no longer cancels the spin- wave contribution. Essentially 



g{r - r') = exp - 



exp 



(^-i3|ln|(r-r>,||'^')|(r-r>p|-V2, 



(80) 



where we systematically dropped all contributions going slower than any power law. We see that, contrary to what 
happens in the usual case, the total correlator is not of the Ornstein-Zernike type. However, in the limit _B — > 0, the 
Ornstein-Zernike form is retrieved. 
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F. Analytic continuation 

It is now necessary to obtain the phase correlator in real time, and then in real frequency. We need to perform an 
analytical continuation t — > i t. We do not attempt to get the short time behavior. This would correspond to high 
voltages. At long times, 

p{t) cx eyi^{-\ujpt\C'^)\ujpt\-^l'^ exp(-B |ln |a;pt| |^^^) , (81) 

where cx means "proportional to". Normalization of p would require more information on the shape oip{t) for shorter 
times. This yields in real frequencies f2, for small 17 's (with respect to ujp) 

P{n) cx e(l]-Wpr^)exp(-B|ln|f]/wp-r^l|^^^) (82) 

where Q is the Heaviside function. For larger f2, this form is not valid. One would need the short time behavior of 
p{t), which we do not have explicitly. However, on those short timescales, the long range part of the interaction will 
play a minor role. Therefore, in order to perform calculations, a reasonable assumption is that, as without LRI, P(r2) 
is negligible^^ for $7 larger than nujp, with n = 2. Enforcing the normalization P{i})dV, — 1, we thus have the 
approximate form 

with po = [2C{2^/'^B, ^/ k — ^^^)] ^, and the function C is defined by C{b, y) — exp( —6 In \x\ ] dx. 



p{n) = pouj^'ein/Ljp^c') e^p[-B\in\n/LOp- r'\\ ' ) {^/l^p- CT'^', (83) 

3/2\ 



V. PHYSICAL APPLICATIONS 

Having obtained the phase correlator in real time, we then use it for the calculation of the Andreev current. A 
normal metal is placed on the left of the Josephson junction chain and biased by a voltage V, below the gap |A|, see 
Fig. 1. 

A. Andreev average current 

The average Andreev current flowing from the metal into the superconductor can be readily expressed as a function 
of the distribution of phase modes P{^)- The case P{^) = S{^) corresponds to a usual BCS superconductor, and 
BTK results^^ are retrieved. The case where P{ft) is given by the Luttinger liquid theory or usual 2d XY model 
has been studied in Refs. HHi. The contact between the metal and the JJ chain can range from tunnel to perfect 
contact. Following the Keldysh formalism, and performing a nonperturbative expansion, to all orders in the tunneling 
amplitude, one can describe all situations'^. Without phase fluctuations, the procedure has been shown to be exactly 
equivalent to solving the Bogolubov-de Gennes (BdG) equations. With phase fluctuations, it is only an approximation. 
The average current-^ reads: 

IfipF^ P ~ /""R — 2" 

^^^ = Tm^^^^ Jo ^ Hrf^, (84) 

with 

(^--^11 1 + ^1 +4^-^ (85) 

W being the bandwidth and F a hopping parameter. For tunnel contacts, |r| ^ W. A perfect contact corresponds 
toT = W exactly. See Ref. [H for details. 

We would like to compare two cases, with and without LRI. For JJ chains, LRI are usually not taken into account. 
One exception are the models of JJ chain considered in Ref. l27i . where long range interactions in the capacitance 
matrix were examined. These models yield physical results which are esentially the same as the usual XY model. 
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except in the case where capacitance of the islands to the ground, is strictly zero. Only in this case does the 
BKT transition disappear and the chain is always in the disordered phase, because of the proliferation of unbound 
vortices. 

In most experimental situations, it does not seem to be realistic to assume that C^^ is really completely negligible 
with respect to nearest-neighbor islands capacitive coupling (or even with respect to longer range inter-island capacitive 
couplings). One exception may be the case of the experiments of Ref. 6. By contrast, the model studied in this paper 
always keeps a non-zero capacitance to the ground and still always gives unbound vortices. 

Now, we would like to compare the results of our calculation with a model commonly used in the literature. 
In this latter model, a nearest-neighbor interaction bewteen island charges is included, and there is no LRI. We 
discuss qualitatively this case and compare the predictions from such a model with our case. A term of the form 
X^i 9i9j+i/(2Ci) must be added to the Hamiltonian, see Eq.®. Qualitatively, this will add some charge fluctuations 
and should decrease the phase correlations. In a more quantitative way, this will change for example the correlator of 
6 in Fourier space Eq. ([T4)) to 



ttK 



■^q,uj^ — q,—uj/ 



q^[l + §^ cos(ga)] 



(86) 



where a is the spacing between Josephson junctions. Calculations are straightforward and the main changes with 
respect to the case with just a local Coulomb term are the following. In the quasi long range phase, for frequencies 

much smaller than cj„ 



1-t- 



and thus the constant K has to be changed 



Up, the plasma frequency is now — w. 

accordingly to K' = K^l + This no longer holds for frequencies of order ujp. The critical point will be shifted 
Defining t = hup/J, the usual BKT transition which occurred at t = tc ^ tt/2 will now take place at t'^ 



In the disordered phase, but still no too far from criticality, the correlation length will be depressed. However, the 
phase correlator will retain approximately its Ornstein-Zernike form, Eq. ()78|) . with different values of Up and f . The 



important thing is that the pre-exponential factor is still in 1/ \/\v — r'|. 

As an example, we take a perfectly transmitting interface and a situation where the Josephson chain would be 
insulating (disordered phase) even in absence of LRI. The presence of LRI will modify the correlation length ^ but 
also alter the form of the correlator P{^). Then, D{ijj) is a constant and the current assumes the form 



{I) 



Po 



X 



exp 



(^-B\ In x\ 



3/2 



{X — x) dx, 



(87) 



with 



X = 



2eV 



2g ^ ^, just as without LRI. The new thing is that, for voltages slightly 



(/) remains zero for voltages smaller than Vq 

larger than Vb, (/) no longer behaves as {V — Vb)'^/^. We introduce the dimensionless quantity = exp 
and the corresponding voltage Vi = -^{l^^ + C^^). Vi is a crossover value of the voltage, not a sharp threshold as 
Vq. For voltages V, such that Vq < F ^ Vi, then, X <^l^^ and (/) behaves as 



3/2 



(89) 



with 6*2 a constant. Thus, the current starts to increase faster than any power law of {V — Vq). If, on the other 
hand, V is larger than Vi, then, (J) increases again as {V — V^f/'^. For even larger voltages, properties of the phase 
correlator at short time scales are needed. We expect that, the LRI will not be relevant on these scales and that 
the linear behavior of the I ~ V curves in the usual case will be retrieved. This is illustrated in Fig. 5, where we 
plot (/) in units of ewp/vr versus eV/{2hujp), for t = 27r. Since our calculations describe only moderate LRI, for the 
LRI parameter A, we take A = 0.2 which results in S = 0.316. For the inverse correlation "length" we take^ 

it to be exp^— 1.5/V^ — ■ This serves just as a rough estimate since the exponential behavior of the correlation 

length applies only (above and) close to the usual BKT transition point. It leads to ^^^/4 ~ 0.125, to be compared 
to /~^/4 ~ 0.029. We see that the Andreev current is strongly depressed by LRI. We also compare with the model 
without LRI but with nearest neighbor Coulomb interactions, mentioned above. We have taken ^ = 0.25. This 
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causes the voltage threshold Vq to be increased to a slightly higher value Vq, but the current still starts as {V — Vq)^^^; 
see Fig. (5). 

To be more quantitative, we make a logarithmic plot, on Fig. 6, of (/) versus [{V — Vo/Vq)]^^^. Whereas the curve 
without LRI, for A ~ has a slope 3/2, this is not the case with LRI. The solid curve catches only with the curve 
without LRI, for y — Vb of the order of 0.1 to 1, which is slightly larger than of order /~^/4. For the model with 
nearest-neighbor interaction, we plot instead (/) versus [{V — Vq)/Vq]^' The curve (dash-dotted line) is still very 
close to a straight line of slope 3/2. 

In the case where the Josephson chain would be in the (quasi) ordered phase, t < tc, without LRI, LRI impose a 
correlation length of order 

^1 = A~^vk^, (90) 

and thus a threshold voltage {hL0p/2e) ^. 




FIG. 5: Andreev current (7), [ in units of CLOp/n ], as a function of applied voltage V, [ in units of 2fkOp/e ], for a perfectly 
transmitting interface and Hup/J = 6.28, with LRI and A = 0.2 (solid line), without LRI {A = 0) but with a nearest-neighbor 
interaction (dash-dotted line), with only on-site interaction (dashed line). 




10"^ io"^[(V-Vo)A/o)]^ 10"' 1. 

FIG. 6: log — log plot of Andreev current {/), as a function of [(V - Vo)/Vo^ ^ ^, (see text), same legend and values as Fig. 5, 
except for the case with a nearest-neigbor interaction (dash-dotted line). In this case the threshold value Vq is different from 
the threshold Vo for the two other cases and we plot {V — Vo)/Vo ^' ~^ instead on abscissa. 
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VI. CONCLUSION 

A study of the influence of unscreened long range Coulomb interaction was carried out for one-dimensional array of 
Josephson junctions. This can also have implications for one-dimensional thin superconducting wires. We focussed on 
the case where the initial term describing LRI was much smaller than Josephson coupling and local Coulomb blockade 
term. Deriving a Hamiltonian describing the interaction of phase slips, we examined how the KT flow is modified by 
LRI. 

Implications on Andreev reflexion between a normal metal and a superconductor was examined and LRI introduce 
a new voltage scale. The analytical form of the I — V characteristics was derived for any value of the transparency of 
the normal-superconductor interface, (from tunnel contact to perfect interface) . If we start from a situation where the 
charging energy dominates the Josephson coupling, corresponding to the JJ chain being insulating, the introduction 
of LRI changes the behaviour above the threshold voltage Vc from a {V — VcY^^ behavior to a sharper increase, faster 
than any power-law of {V — Vc). This lasts up to a crossover voltage Vi, whose value depends on the strength of 
the long range Coulomb interactions. If, on the contrary, we start from a situation where the Josephson coupling 
dominates without LRI (quasi ordered phase), then LRI will always generate a threshold voltage Vq below which no 

Andreev current can flow. This threshold Vq goes as A vt^-t where A is proportional to the strength of the LRI and 
t — tc is the distance to the critical point in the absence of LRI, with t — hujp/J and = tt/2. With LRI, above 
Vq, the current starts to increase faster than any power law. The form of the phase correlator obtained in this paper 
could be used for calculation of other physical quantities relevant to JJ chains or to one-dimensional superconducting 
thin quantum wires. 

Such calculations are not directly relevant for high- Tc superconductors^, because of their layered structure and 
their strong anisotropy in the direction transverse to the planes, whereas the former calculation remained striclty 
one-dimensional. 

For designed JJ chains, the inverse capacitance is generally cut exponentially with a screening length of the order 
of As ~ aJC jCg., with usually As ~ 30a, (a is the spacing between neighboring junctions),— i^. The number TV of 
junctions^ can be up to 255, so that Xs <ti L = N a. Thus, the model considered here should not be construed as a 
model for fabricated JJ chains but rather an attempt to include LRI, which do exist in isolated thin wires. However, 
we can use the above calculations, replacing the badly screened Coulomb potential by the exponentially screened 
potential Uij , see Eq. ^ . This leads qualitatively to an interaction between vortices Vp p, which has the following 

properties. For distances p much smaller than Aa, Vp,p' = 2?y^^^A^^/^(a//9) Inp, with T) a constant of the order of 

ten. The strength of the interaction between vortices decays very rapidly, (much faster than for the badly screened 
Coulomb interaction). This lasts till p becomes of the order of Aa, where V'p p, stops to decrease and behaves as 

Vp^p' — JE^^ '^^'^ Lnp, for large distances, i.e. p >• aA, 

Theoretically, the RG group is started at coarse graining scales smaller than Aa and as the RG flow proceeds, 
the interaction between vortices gets smaller. This may tend to drive towards the disordered region. This is not 
unexpected, given the large amount of Coulomb repulsion provided by the charge coupling of the islands which are 
less than Aa apart. When RG proceeds further and the coarse-graining length becomes larger than Aa, then, the 
interaction Vp^pi bewteen vortices ceases to decay and becomes logarithmic with distance. According to the value of 
the coupling reached at that moment, we can flow either to the disorded region or to a quasi-ordered flxed point. The 
screening of the Coulomb interaction renders quasi long range-ordering more favorable. In particular, contrary to the 
badly screened interaction case, phase correlations do not always decay exponentially at the end. 

Let us suppose first that J is sufiiciently large so that we end in the (quasi)-ordered phase, P(il) will decay as a 
power law for small fl, with some exponent a, for f2 smaller than uJp/X. For Q larger than ujp/X, the decay of P{^) 
corresponds to a larger effective vortex interaction, and phase correlations are enhanced, leading to a larger P{^). 
As we go to large distances (or longer times or smaller voltages), vortex interactions get reduced until the distance 
Aa is reached. This is due large Coulomb repulsion at shorter lengthscales. This has experimental consequences on 
the shape of / vs. V curves. Frequencies smaller than ujp/X correspond to voltages smaller than hiUp/e. For these 

voltages, I — V curves will start as power law with some exponent, normally as V" with a = 1 + ■^), where 
Eq and J' are renormalized values of Eq and J. When V becomes larger than fuop/e, the current is enhanced and 
increases faster than V"' . This remains true for voltages still smaller than hujpj e. 

Now, let us suppose that J is not so large and one ends in the disordered phase. There is a threshold voltage but 
beyond Vc, the current starts to increase as {V — Vcf'^'^ in the usual theory. Inclusion of screened interaction leads to 
an increase which is faster than {V — Vc)^^^, which could also be tested experimentally. 

The model with badly screened interactions may also have applications for very long (but not unrealistically long) 
Id superconductor thin wires. More precisely in thin wires (typically 3 to 5 nm diameter), the interaction between 
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phase slips has the same structure as the Hamiltonian (1), but the physical origin is somewhat different—. First, in 
concrete situations, it is not the two-dimensional BTK transition which is seen but rather a one-dimensional Schmid 
transitional. We briefly recall the reasons below. 

In experiments with very homogenous thin wires in the difusive limit, of diameter 3 to 10 nm, such as the ones used 
in Refs. 29,30 the maximum length L is about 1/im. Operating temperatures are between lOOmK and a few Kelvin. If 
the length were infinite, renormalization had to be stopped when coarse graining length r on the imaginary time axis 
reaches I3h. However, in practise, two-dimensional renormalization has to be stopped when it reaches L/vms where 
VMS is the Mooj-Schon velocity^, vms ~ ujpd/CL is equivalent of the spin- wave velocity, where Cl is the capacity per 
unit length and d the diameter of the wire, vms is of the order lO^m.s"^, so that for a 1/xm long wire, L/vg is of the 
order 10~^^s and for temperatures lower than 0.75K, renormalization has to be stopped first on the x axis, when the 
coarse graining length reaches the size of the wire. We end up with charges in 1 dimension (as for a single Josephson 
junction) which interact logarithmically as IrenLn\ri — rj\, with Iren a renormalized coupling. Then, renormalization 
proceeds solely on the imaginary time direction. Depending on the value of Xren , fugacity of vortices renormalizes to 
zero or not, resulting in the usual Schmid transition. The spatial structure on the x axis has been completely lost and 
the long-range interaction do not play any role, except that they changed slightly the values of I, en. Thus, in this 
case, LRI do not seem to play an important role. Thus, in order to see the influence of LRI, longer wires are needed. 

For longer wires (but not totally unrealistic, a few microns long), we can expect temperatures larger than hvsl {Lks), 
so that renormalization has to be cut on the imaginary time axis flrst. Then, LRI will drive the resistance i? to a 
flnite value. If R versus T curves arc plotted, no curve should plunge towards zero resistance, unless T becomes of 
the order of hvs/ (LkB), which is about 0.5 K for L = lO^m. 

There are difficulties to see LRI. One could be that stray capacities (unavoidable in experiments) can also cut the 
long-range term in the capacitance for small q (in real space) so that C^^{k) is no longer of the form Infc for small 
k. This depends on the charateristics of the wires. Another problem is that the action for a single phase slip, which 
controls its bare fugacity, is strongly dependent on temperature and on the diameter of the wire. This dependence 
seems to be more important than the interaction bewteen phase-slips in determining the shape of the curves showing 
the resistance R versus T, for different wire thicknessesSS. 

Long range interactions are originally not screened in isolated one-dimcnsional wires. However, the presence of an 
environment could alter the screening. Situations where a plate or a wire is placed in the vicinity of a superconducting 
thin wire are now experimentally?^^ and theoretically studied^i^. This aims at providing mechanisms of dissipation 
so as to damp the phase fluctuations. Though we did not consider in this paper coupling to an external environment, 
our calculation may be viewed as a first step towards understanding the role of LRI, including, both spin waves and 
phase slips (solitons). 
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Appendix A 



In this appendix, we derive the RG flow in the presence of long-range interactions. We define g as 

g{R,,,) = JJ^{r,~rj), (Al) 

with 

R,., = \v,-r,\. (A2) 

Integration has to be done over in the whole plane except the two circles of radii t around and r;. Setting 
r = (r; + rfe)/2 and denoting by /C the vortex fugacity, the partition function Z reads 



T^(rj~rfc)(rj-r;) 
\rj - rfc|2 \rj -r,|2 



exp|^ J'(ri -rj)gjgj ln|r, -rj|,|. 



(A3) 
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where D is the whole plane except aera inside the circle of radii r around the r^, i = 1, n. Let us denote by 

where T) is the whole plane except the two circles of radii r around rfc and r;. Given the value of J-'ir), see Eq. (|40p . 
we have 



(A4) 



I ~ 27r HRk,i/T) [1 - -A^^HRk,i/T) 



for <C 1. Denote by 



C = 27rT dr 



k,l k^l 

where S is the aera. Collecting all the terms, 
Using = 1 + x, for small x, 



S - 27rT2^2|^ _ -^0n(i?fc,,/r)}gfc9, In 



rfc - I"; 



n+l ri+1 



^+ (^ + 1)2 exp(^^5(i?,,,)q,;g, 



i=l j=l 



In 



r,: - r,- 



n+l n+l 

Z => ^/C-|exp{^^,_^^, 



(n + l)' 



j (fri... J (fvn exp J ^ g{Ri,j)qiqj in 



27rTdT5 



(A5) 



(A6) 



(A7) 



dr 



.(2^)V^/C^ — ^^^^ 

k.lk^l 



R 



k.l 



1 - -A\IHRuj/t) 



(A8) 



When changing r into r + dr, this leads only to the renormalization of the fugacity JC into K,{l — J replacing 

IF{ri — Vj) by its initial value, which is one. 

After changing (fc, /) — (*, j) in the second summation, the part in the exponential reads 



(A9) 



If J-{Rij) were a constant, we would retrieve the usual KT flow. However, since we are looking at larger and larger 
lengthscales, the effective coupling will not only decrease because of RG but also because of its Vln R dependence. 
The smallest R contributing to Z is r. Thus, adding a term proportional to —^A\Jh\[T) will take this into account. 

Now, the partition function has to be rearranged so as to look as its original form, with renormalized parameters 
and an overall multiplicative factor, setting z = In r, we obtain 



dj_ 

dz 
die 
dz 



-JJC. 



Switching to the usual variables 



and linearizing around A" = 0, the flow reads 



X 
Y 



dX 

dz 
dY 

dz 



^-2, 

-Y - BV^, 
-2XY, 



(AlO) 
(All) 

(A12) 

(A13) 
(A14) 

(A15) 
(A16) 



with B = 4\/A/3. 
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Appendix B 

Here, we give details on the contributions to the vortex part of the correlation function ^^(r — r') coming from 
smooth variations of v{r). The first term to evaluate is 

T[ = - j j \v{K)(fn, (Bl) 

where the double integral extends over the whole shaded surfaee in Fig. 4. A is the Laplacian operator (with respect 
to variable R). In this calculation, r' and r are fixed, and, without restriction, r is taken to be zero for simplicity. If 
R is the modulus of R, since v decreases as ln(i?)~^/^ for large R, it turns out that Av behaves, as R~'^ (In J?)~^/^. 
Thus, V /S.V behaves as (lni?)~^. When surface integration is performed, the integral tends to a constant, for 
r' — > oo because of the (lnii)~^ factor. When exponentiated this will give a constant. Since, we shall neglect all 
terms that depend on r', weaker than any power law, we shall neglect this factor. 
The second term is a line integral over the outer contour Tout on Fig. 4, 

^ I v{\/v).da, (B2) 

where da is the normal to the contour Tout which is here chosen as a circle of radius TZ. For all R's belonging to Tout, 
R are equal to TZ and v are of order (In 7^)~^/^, the leading terms in iVv) are of order TZ~^{\a TV)~^^'^ . Thus T2 is at 
most of order (In TZ)~^ and vanishes as TZ goes to infinity. 

The third term to evaluate is a line integral over the inner contour Fj^. on Fig. 4. 

T^= [ v{Vv).da. (B3) 

A calculation completely analogous to the case without LRl shows that the leading terms are of order 
/:' ^{ln(.- + y-)}-' d.' -H ^ (.' - with r' the modulus of r - r'. This behaves, for large r' as 
ln(ln r'). 
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